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ABSTRACT 

We explore whether protoplanetary disks with self-shadowing from puffed up 
inner rims exhibit observable features in scattered light images. We use both 
self-consistent hydrostatic equilibrium calculations and parameterized models to 
produce the vertically puffed up inner rims. We find that, in general, the transi¬ 
tion between the shadowed and flared regions occurs in a smooth manner over a 
broad radius range, and no sudden jump exists at the outer edge of the shadow 
in either the disk temperature or density structures. As a result, a puffed up rim 
cannot create sharp ring/arc/spiral-arm-like features in the outer disk as have 
been detected in recent direct NIR imaging of disks. On the other hand, if the 
puffed up rim has a sharp edge in the vertical direction, the shadowing effect can 
produce a distinct 3-stage broken power law in the radial intensity profile of the 
scattered light, with 2 steep surface brightness radial profiles in the inner and 
outer disk joined by a shallow transition region around the shadow edge. These 
types of scattered light profiles may have already been observed, such as in the 
recent Subaru direct imaging of the TW Hydrae system. 

Subject headings: protoplanetary disks — stars: pre-main sequence — stars: 
variables: T Tauri, Herbig Ae/Be — circumstellar matter — radiative transfer 
— planets and satellites: formation 


1. Introduction 


Flattened, rotating circumstellar disks reprocess light from newly born stars in T Tauri 


and Herbig Ae/Be systems (Appenzeller & Mundt 1989; Waters & Waelkens 1998). Disks 
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can modify the spectral energy distribution (SED) of the systems, and reveal themselves 
in resolved images at various wavelengths (Williams & Cieza 2011). In a disk, dust grains 
only exist outwards the sublimation radius r su b, where dust reach sublimation temperature 
T sub , leaving a dust-free hole at the center. Protoplanetary disks are usually optically thick 


to stellar radiation (with exceptions such as transitional disks, Espaillat et al. 2014). As 
a result, most part of the disk receives direct starlight only at the surface, under a small 


grazing angle determined by the shape of the surface (Chiang & Goldreich 1997; Takami 


et al. 2014). On the other hand, the entire wall of the rim at r sub is in full sight of the 


star, heated by intense stellar radiation. As a result, the rim is vertically puffed up under 
the hydrostatic equilibrium condition as pressure balances gravity (Isella & Natta 2005). As 


shown by Dullemond et al. (2001) and Dullemond & Dominik (2004a), a puffed inner rim 


blocks starlight from reaching to the disk behind, thereby creating a shadow in the outer 
disk. Depending on the disk properties, the outer disk may re-emerge from the shadow at a 
certain distance, or remain shadowed throughout. The puffed up rim and the self-shadowing 
effect are more prominent in Herbig Ae/Be stars than in T Tauri stars, as dust sublimation 
occurs closer to the less luminous and cooler T Tauri stars, and the covering angle of the 


rim is smaller (Dullemond et al. 2001). 


Self-shadowing by a disk rim is a key component in explaining observations of Herbig 


Ae/Be stars. In particular, the difference between group I and II objects (Meeus et al. 
2001]) has been interpreted as that group I disks are flared while group II disks are flat and 


(completely) self-shadowed by a puffed up rim (Dullemond & Dominik 2004a; Dullemond 
et al. 2001]). Also, a puffed up rim has been involved to explain the light curves of UX Orionis 
stars (Dullemond et al.||2003), dust emission features in disks (Acke fc van den Ancker|2004), 


and cooling of the surface layer in disks (Brown et al. 2012). The change of rim height has 


been proposed to account for the variability of disk infrared excess (Espaillat et al. 2011 


Muzerolle et al. 2009; 

Juhasz et al. 2007 

) as well as the overall surface brightness of the disk 

in scattered light imaging ( 

Wisniewski et al. 2008 

). Moreover, an evolutionary sequence has 

been proposed by Acke et al. (2004) and 

Dullemond & Dominik 

(2004a): as the shadow cast 


by the rim extends further as the amount of small grains decreases due to grain growth and 
settling, eventually an initially flared group I object evolves into a flat self-shadowed group 
II object at late stages. 

Recently, thanks to a few high angular resolution imaging instruments that have been 
mounted on 8-meter class ground based mirrors in the past few years, near-infrared (NIR) 
imaging of polarized scattered light from disks have found structures that may be produced 
by self-shadowing in disks. Marino et akl (2015) suggested that the two nulls seen on the 


ring in the H band polarized intensity (PI) image of the HD 142527 system are shadows cast 
by an misaligned inner disk (see also |Verhoeff et al? 2011, Avenhaus et al. 2014, Fukagawa 











































































et al. 2013, Perez et al. 2014, and Casassus et al. 2013). Grady et al. (2013) found that 


the two spiral arms in the MWC 758 system may have shadows in its H band PI images 
(see also Garufi et al. 2014). In addition, various rings and arcs and gaps in disks have also 


been suggested to be manifestations of self-shadowing (e.g. HD 144432, Chen et al. 2012 


HD 169142, Quanz et al. 2013 ; HD 163296, [Garufi et al.| 2014 ; TW Hydrae, |Akiyama et al. 


2015; HD 169142, Momose et al. 2015). In these cases, the gap region has been interpreted 


as being shadowed by puffed up inner disk structures, with the outer gap edge presumed to 
be the edge of the shadow. 


On the theory side, the effect of shadowing have been explored by many in the past 
through radiative transfer modeling, with the shadow cast by either a puffed up inner rim or 


hydro features produced by planets (e.g. 

Dullemond et al. 2001; 

Dullemond 2002 

; Dullemond 

& Dominik 2004a[ van Boekel et al.|2005| Meijer et al. 2008| Tannirkulam et al. 2008[ Jang- 

Condell 

2009; 

Jang-Condell & Turner 

2012; 

Siebenmorgen & Heymann 

2012) 

. However, 


most works focused on the SED or the thermal emission of the disk. The effects of a puffed 
up rim in high angular resolution scattered light imaging have not been throughly examined. 
In particular, whether observed rings and arcs in disks can be the outer edge of shadows 
remain unanswered. 


In this paper, the properties of the shadow cast by a puffed up inner rim in NIR polarized 
light imaging are studied through radiative transfer simulations using the Whitney et al. 


(2013) code. Disk vertical structures are self-consistently calculated assuming hydro static 
equilibrium (HSEQ). Synthetic images of disks with either a naturally puffed up rim or 
a parametrized rim are produced, and their properties are measured and compared with 
observations. The structure of the paper is as following. The models and MCRT simulations 
are introduced in Section [2| the results are presented in Section [3j followed by a summary 
and discussion in Section 01 


2. Simulation Setup 


MCRT simulations are carried out using the Whitney et al. (2013, see also Whitney 


et al. 

2003a|b 

code, which has been used to model protoplanetary disks ir 

i the past (e.g. 

Hashimoto et al. 2012; Zhu et al. 2012; Dong et al. 2012b ; a[ Follette et al. 2013; Grady 

et al. 2013 

. The simulation setup is largely adopted from 

Dong et al. 

(2012b 

), and is briefly 


summarized here. The central star has M* = 2.5M 0 , R k = 2R & , and T k = 10 4 K, typical for a 
Herbig Ae/Be system. The disk is axisymmetric, and is constructed in spherical coordinates 
r (radial) and 6 (polar, midplane is at 0 = 0 and poles are at 0 = ±90°). The grid has 
500 cells in r direction with logarithmic spacing and 200 cells in 0 direction with power law 
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spacing to over sample the inner disk and the disk midplane. The inner edge of the disk is 
at dust sublimation radius r su b ~0.5 AU, where dust temperature reaches 1600 K. The outer 
edge of the disk is at 200 AU. The surface density S follows a power law with radius as 


£(r) = E, 


r 0 


( 1 ) 


where E 0 is the normalization at r 0 . The 1/r power law has been suggested by both detailed 


disk structure calculation (e.g. D’Alessio et al. 1998) and mm disk observations (e.g. Andrews 


et al. 2009, 2010). The total mass of the disk is assumed to be = 0.01 M* = 0.025M o . 
A small accretion rate of lO _9 M 0 /yr is assumed to represent a typical accretion state, and 
also to provide nominal heating to the shadowed region to prevent the development of a 
known ripple instability (Dullemond 2000; [Chiang 2000). All simulations are run with 1 
billion photon packages. Synthetic H band (1.6 pm) image^j] and SED are produced for 
each model. All models are assumed to be at 140 pc away. Full-resolution MCRT images 
are convolved by a circular Gaussian kernel with a full width half max (FWHM) of 0.05" to 
achieve an angular resolution comparable to observations by Subaru, VLT, and Gemini (the 
FWHM of an airy disk is 1.028A /D ~ 0.04" for a primary mirror with a diameter h = 8.2m 
at A = 1.6 p m). 

In simulations, photons from the central star are absorbed and reemitted or scattered by 
the dust in the surrounding disk. The temperature in each grid cell is calculated based on the 


radiative equilibrium algorithm described in Lucy (1999). The anisotropic scattering phase 


function is approximated using the Henyey-Greenstein function. Polarization is calculated 


assuming a Rayleigh-like phase function for the linear polarization (White 1979). The vertical 


density structure of the disk is solved through iterations assuming self-consistent HSEQ in 
the vertical direction z, as gas pressure balances gravity: 


dp 

dz 


zGM * 
P 2 ' 


( 2 ) 


where p is pressure and p is density. Simulations are iterated on the temperature T calculated 
from the radiative equilibrium solution of the radiation transfer equation at least 10 times 


to achieve convergence on the vertical structures (Whitney et al. 2013). In this study the 


stellar radiation is the main energy source of the disk (heating from accretion has little effect 
except in the umbra region right behind the rim), and the gravity from the central star is 
the sole gravitational source. Note that the vertical structure of the disk (as well as the 


1 The physical quantity recorded in all images is the specific intensity, or intensity I v for short, which 
has the unit [mJy arcsec -2 ], or [ergs s -1 cm -2 Hz -1 arcsec -2 ]. This quantity is sometimes referred to as 
“spectral radiance”. 
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puffed up rim) may be affected by other mechanisms, such as the gravitational perturbation 
from an orbiting companion (particularly if on an inclined orbit), or turbulences in the disk. 
The naturally puffed up rims based on the HSEQ condition in this study do not account for 
these additional mechanisms. 


There are two types of dust grains in the models, the “small” grains and the “big” 
grains. The gas-to-dust-mass ratio (including both the small and big grains) is assumed to 
be 100:1 (therefore the total mass of the dust grains m gra i ns is 2.5 x 1O _4 M 0 ). The small 
grains represent primordial dust particles. Their properties are assumed to be the standard 
interstellar medium (ISM) grains as in Kim et al. (1994). These grains contain silicate, 


graphite, and amorphous carbon, and their size distribution is a smooth power law in the 
range of 0.02-0.25 /mi followed by a sharp cut off beyond 0.25 //in. The small grains are 
well mixed with gas as the two are dynamically well couplec0 The large grains represent an 
advanced stage of grain growth and coagulations. Model 1 in Wood et al. (2002) is adopted 
for their properties. They are made of amorphous carbon and astronomical silicates, and 
have a power-law size distribution as n(s) ~ s -3 , an exponential cutoff beyond 50/im, and 
a maxim size of 1 mm. The big grains tend to settle to the disk midplane, as they usually 
dynamically decouple from the gas and are no longer held in the vertical direction by the 


gas pressure (Dullemond & Dominik 2004b). The vertical distribution of the big grains 
is assumed to be a Gaussian profile with a small scale height hb g = 0.01r(r/lAU)°' 25 (i.e. 
their vertical structure is parametrized and NOT evolved in the HSEQ iterations). Therefor 
the big grains form a thin layer at the disk midplane. The details of the properties and 
spatial distribution of the big grains have little effect on the scattered light images and 
SED presented below, as long as they are well embedded at the disk midplane. The optical 


properties of both grains can be found in Figure 2 in Dong et al. (2012b). 


Two sets of models, series A and B, are constructed and listed in Table [lj Model Base 
is the fiducial model, which has no big grains. Series A is constructed by fixing m gra i ns while 
decreasing m gra i nS)Sma n and increasing m gra j nS; bi g , similar to the ML series in Meijer et al. 


( |2008| ) and BL series in Dullemond & Dominik (2004a). This scenario (from model Base to 
A4) is to mimic an evolutionary sequence: the amount of small grains decreases with time 
as grain growth and coagulation gradually convert them into bigger particles that settles 
down at the disk mid-plane. The inner rim in Base to A4 is naturally puffed up under 
the HSEQ condition (only the small grains). As will be shown in the next section, while 


2 It is possible that /jm-sized grains may still have some degree of settling, depending on the strength of 
the turbulences in disks that stir up dust grains vertically (Dullemond & Dominik 2004b). Due to a lack of 
good understanding of the nature and strength of the turbulences in protoplanetary disk, settling of these 
sul:>-//in-sized small grains is not considered in this work. 
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series A does gradually turn from a flared disk (Base) into a self-shadowed disk (A4), the 
boundary between the shadowed inner disk and the flared outer disk is not sharp. This is 
because the rim does not have a sharp edge in the vertical density distribution, and the radial 
optical depth of the rim only gradually decreases with height. To maximize the sharpness 
of the rim edge and produce features with the highest contrast possible in scattered light 
imaging, series B is created. In B1 to B6, an artificial opaque rim is placed at r = r su b to 
1.1 r su b, with a parametrized height H r i m ,opaque expressed in the angle subtended by hah of 
the rim to the star (i.e. the vertical structure of the artificial rim is fixed and NOT evolved 
under the HSEQ condition). The opaque rim is constructed by filling the region with small 
grains of a uniform density 10 -9 g cm -3 , so that the radial optical depth r of the rim is 
greater than 10 4 at wavelengths of interests, therefore completely blocks the starlight in the 
solid angle subtended by the rim. In other words, in the polar direction dp/d,6 = 0 inside 
the rim, then p = 0 outside the rim. This artificial vertically constant density structure is 
only designed in order to create the sharpest possible shadows. It is not meant to have a 
physical motivation. Both the naturally puffed up rim and the parametrized artificial rim 
are illustrated in Figure [TJ 


3. Results 


This section presents the main modeling results, first the structure of the disks (Sec¬ 
tion 3.1) and then observational aspects of the models (Section |3.2|). 


3.1. Disk Structure 

The density and temperature structure for models Base and B2 are shown in Figure [l] in 
spherical coordinates, with the disk surface over-plotted, defined as where the radial optical 
depth r at 0.4 //in at each polar angle 0 from the stellar surface (the stellar flux peaks 
aroundt 0.4 pm) reaches 1. In spherical coordinates a shadowed region has a flat surface 
(constant 6 with increasing r), while a flared region has a rising surface (increasing 6 with 
increasing r). The disk surface traces out the boundary between the hot disk atmosphere 
and the cool interior, as most starlight intercepted by the disk is deposited near the surface. 

In the model Base, a puffed up inner rim is clearly visible in the density structure as 
a vertically extended structure close to r su | ;) . The rim casts a shadow behind from r su b to 
~ 3 AU, where the surface comes out of the shadow and starts rising. However, as discussed 
in Section [2j the puffed up rim has a blurry edge. As a result, the edge of the shadow is not 
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sharp, and the disk surface smoothly transits from flat to flared around 3 AU. 

In the model B2, a parametrized artificial rim is placed at r su b with H T i m ,opaque = 12.5°. 
The sharp rim edge results in a sharp turning point on the r = 1 surface at r s h a dow = 13 AU, 
where the surface suddenly rises. Inside the shadow region, the disk does not have a uniform 
structure, and can be separated into two regions. An umbra exists immediately behind the 
rim (boundary indicated by the dashed line), within which the entire surface of the star is 
not visible. In the umbra, practically the only heating sources are radiation diffusion and 
accretion. This region becomes sufficiently cold to collapse to the midplane. On the other 
hand, the region in between the umbra and the r = 1 surface receives additional heating 
from the following ways: 


1. Once outside the umbra and moving towards r s h a dow> the disk receives direct radiation 
from a bigger and bigger polar region on the star. 

2. Inward heat diffusion from r s h a dow heats up the adjacent region, smoothing out the 
transition and preventing a sudden jump in T m idpiane (Dullemond et al. 2001). 

3. Indirect illumination and heating from the halo ( |Dullemond fe Dominik 2004a). In 
particular, closer to r s } ia d ow • secondary (multiple) scattering from the disk atmosphere 
in the flared region becomes more important. 

4. NIR flux from the opposite side of the inner rim. 


The strength of additional heating sources 1-3 increases with radius, puffs up the disk mod¬ 
erately starting at ~2 AU. Eventually the r = 1 surface comes out of the shadow at 13 AU 
and disk stays flared outward. 


The t — 1 surface for all models are plotted in Figure [2} As shown in Dullemond & 


Dominik (2004a) and|Meijer et al. (2008), from Base to A4 the disk gradually evolves from 


flared to flat (completely self-shadowed by its own puffed up rim), as the amount of small 
grains decreases. The shadowed-to-flared-disk transition occurs more and more smoothly, 
and v shadow cannot be clearly defined. On the other hand, the surface transition is sharp in 
the B series, as the surface suddenly turns upward at [r s h a dow,#rim,opaque] (marked by dots in 
the right panel and labeled in Table [T| The entire disk in B6 is inside the shadow so r s h a dow is 
not defined.). The transition point rises and moves outward as iTrim,opaque increases, and the 
surface always converges to model Base once out of the shadow, indicating that the flared 
part of the disk in all models are similar. 


The midplane temperature of the disk T m id P iane is shown in Figure [3] for the B series. 
There are two bumps on the curves. The first bump occurs when the disk comes out of the 





















umbra, and the second bump occurs at the shadow edge r s h a dow- In the umbra all models 
have the same Tmidpiane profile (best traced out by model B6). The location of the first 
bump moves outward as fAim,opaque increases and the umbra extends. After the first bump 
the shadow reduces the midplane temperature, consistent with Siebenmorgen & Heymann 


(2012), and models converge again to another profile in cases that the two bumps are not 
too close to each other (B3 to B6). Finally, all models converge to the Base model at 
r > r s hadow> echoing the convergence of the r = 1 surface in Figure [2| and demonstrating the 
similarity of the models in the flared region. The r = 1 surface officially comes out of the 
shadow at r shadow and starts to be directly irradiated by the bulk part of the stellar surface 
(|latitude| > iZ r i m ,opaque)- However, the transition to the “flared profile” occurs much earlier 
than r s hadow • For example, B4 starts to deviate from its “shadow profile” (traced out by B6) 
at ~10 AU, and gradually converges towards its “flared profile” (traced out by Base) until 
the eventual merger at r s h a dow = 52 AU. This echoes the smooth transition in temperature 
and density structures seen in Figure [lj 


Vertical density structure of the disk is determined by the temperature structure under 
the HSEQ condition. Therefore, at radius of a fraction of r s h a dow) the disk already starts 
to deviate from its thin “shadow state” and become modestly puffed up, until eventually 
smoothly joins the flared disk at r s h a dow- This is illustrated in Figure [I] for B2, in which case 
the thin “shadow state” ends at ~2 AU, and in between 2 AU and r s h a dow = 13 AU the disk 
is modestly puffed up, and beyond 13 AU the disk becomes fully flared. 


In summary, the disk in Series A evolves from a mostly flared disk to a self-shadowed 
flat disk, but the shadowed-to-flared-disk transition is smooth, due to the lack of a sharp 
rim edge. In the extreme cases, with parametrized sharp rims in series B, the transition 
on the r = 1 surface becomes sharp, but the transition in the underlying disk density and 
temperature structure is still smooth and occurs over a radius range. 


3.2. Scattered Light Images and SED 


SED for all models are shown in Figure |4} Series A mimics the evolutionary sequence 


from group I to II, as shown in Dullemond & Dominik (2004a). The 1-100 //in infrared 


excess drops as the total amount of small grains decreases and the disk becomes more and 
more self-shadowed, meanwhile the far infrared and mm emissions beyond 100 /im rise due 
to more big grains. For series B, the NIR excess from 1-10 /im comes from the artificial rim, 
which increases as the rim grows higher from B1 to B6. Meanwhile the 10-100 /an excess 
decreases as the shadow edge expands and a bigger disk region falls into the shadow. This 
is similar to the “seesaw” phenomenon in Espaillat et al. (2011), which has been explained 
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by puffed up rims with varying heights. Beyond 100 //m the emissions mainly come from 
the outer part of the disk, which remains flared and has a similar density and temperature 
structure in most models thus the emission is little changed. 

Scattered light images trace structures at the disk surface, set by the dust distribution 
and intrinsically determined by the temperature structure under the HSEQ condition. While 
formally inside r s h a dow> the shadowed region may still receives and scatters NIR photons from 
(1) the polar region of the star (in the region immediately below the radial r = 1 surface); (2) 
the bulk part of the star (high in the disk atmosphere); (3) the opposite side of the rim; and 
(4) secondary scattering photons from the disk halo (in particular at radius close to r s h a dow)- 
If the rim is naturally puffed up, stellar photons may also be scattered by the tenuous rim 
edge into the shadowed region. Therefore, it is expected that the shadowed region to be not 
completely dark, and the surface brightness should converge to the flared disk when getting 
closer to r shadow . 

H band polarized intensity images for model B2 and B4 are shown in Figure [5} In B2, 
a marginal ring/gap structure (i.e. a region with PI— r >0 ) at r s h a dow can be seen at face-on 
viewing angle, inside which the shadowed region has low surface brightness. If we define 
“ring/gap contrast” as the ratio of the azimuthally averaged surface brightness at the peak 
of the ring to the bottom of the gap, the contrast is 1.33 for the raw image and 1.03 for 
the convolved image. At a 45° inclination angle, the ring becomes elliptical, and due to the 
forward/backward scattering effect the front side is brighter than the back side 
et al. 2012). In B4, the shadow edge is not distinctive in all panels. At face-on viewing angle, 
the azimuthally averaged surface brightness decreases with radius throughout the disk (see 
Figure |6j). 

The radial profile of face on H band PI images for all models are shown in Figure [6) 
For the A models, converting small grains into vertically settled big grains gradually turns 
a bright flared disk into a faint flat one. The brightness of the disk drops by 1-2 orders 
of magnitude. However this process happens in a rather smooth manner that the surface 
brightness of the entire disk drops all together, while the outer disk drops faster than the inner 
region. The dependence of PI on r changes from PI— r -1 - 8 in the Base model to PI— r ~ 2 - 5 in 
A4. The surface brightness distribution remains smooth and relatively featureless throughout 
in all cases. In particular, no ring/gap structure exists at the edge of the shadowed region, 
in either the original images or the 1/r 2 scaled images. 

For the B models, the shadowed-to-flared-disk transition is clearly visible on the radial 
profiles. In general, the surface brightness radial profile can be described by a 3-stage 
broken power law: (1) well inside the shadow different models have a similar radial profile as 
PI- r~ 1A (best seen in B3 to B6 as their shadows extend far enough); (2) in the flared region 


(Hashimoto 
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all models converge to PI~ r _L8 , the same as the Base model; and (3) a transition region 
with PI ~ r° (constant PI) connects the two. This is consistent with the trend seen earlier 
in the disk density, temperature, and surface structures, that B models have essentially the 
same structure both in the inner shadowed region and the outer flared region, with a wide 
transition region in between. And similar to the T mif i p i. mc profile, the transition occurs over a 
radius range and not suddenly. For B2 to Bfj^J radial profiles smoothly deviate from the r _L4 
profile around r s h a dow/2, and gradually converge to the ~ r -1 - 8 profile around r s h a dow- B3 
has a marginal “gap/ring” structure similar to B2, with very weak ring/gap contrast ratios 
(listed in Table [I]). PI never increases with radius throughout the disk in B4 to B6, and a 
ring/gap structure is not present around the shadow edge. The difference between B2/B3 
and B4/B5 is probably caused by the fact that the transition region is wider in B4/B5 due 
to their large shadow size, so the transition occurs even more smoothly. 


4. Summary and Discussion 

The effects of puffed up inner rims in protoplanetary disks in both temperature and 
density structure and NIR scattered light imaging are studied. Two types of rims are ex¬ 
plored: natural rims that are self-consistently puffed up under the hydrostatic equilibrium 
condition (series A), and artificial opaque rims with parametrized heights and sharp edges 
in the vertical direction (series B). The main results from the models are: 


1. An inner rim naturally puffed up under the HSEQ condition creates a shadow behind, 
within which the disk surface is flat and beyond which it is flared (the A models). De¬ 
creasing the amount of small grains gradually turn the entire disk into a self-shadowed 
and flat disk, confirming previous results such as in |Dullemond et a~ (2001); Dullemond 
& Dominik (2004a); Meijer et al. ( 2008| ) (Figure [2]). 


2. For naturally puffed up rims, the shadowed-to-flared-disk transition occurs smoothly, 
and no sharp/distinct turning point on the disk surface exists. This is mainly due to 
the fact that HSEQ rims do not have a sharp vertical edge, and their “radial optical 
depth” gently decreases with height. As a result, the scattered light images have 
a smooth/featureless surface brightness distribution. From Base to A4, the surface 
brightness of the disk drops by 1-2 orders of magnitude, and the radial dependence 
changes from a shallower r _L8 to a steeper r -2 - 5 (Figure [6]). 


3 r s hadow in B1 is too close to the star and beyond the capability of current NIR imaging facilities. 
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3. The most extreme rims, which are parametrized artificial rims with a razor sharp 
vertical edge in the B series, do create a sharp turning point on the radial r = 1 surface 
(Figure [2]). This point formally defines the edge of the shadow r shadow, beyond which 
the surface is flared and sees the bulk part of the star. Both the inner shadowed region 
and the flared region are similar in models with different rim heights. On the other 
hand, the shadowed-to-flared-disk transition does not occur abruptly in the underlying 
disk density and temperature structures. The midplane temperature starts to deviate 
from its “shadow profile” at a fraction of r s h a dow (Figure [3j), and eventually converges to 
its “flared profile” around r shadow- The density structure of the disk follows the change 
in T m idpi a ne under the HSEQ condition, thereby exhibits a broad transition region as 
well (Figure [I]). This is mainly caused by partial illumination by the polar region of 
the star, multiple scattering, and radiation diffusion within the disk. 

4. Even the parametrized razor sharp rims in the B models cannot produce prominent 
rings and arcs in scattered light images as what have been directly imaged recently. 
On the other hand, self-shadowing by the rim does have a clear effect on the radial 
profile of the scattered light. The surface brightness radial profiles of the B models can 
be described by a 3-stage broken power law: two steep profiles in the inner shadowed 
region and the flared region, connected by a shallow profile in a broad transition region 
(Figure [6]). The enhancement around the shadow edge however is usually not strong 
enough to make a prominent ring/gap structure. B2 and B3 show a very marginal 
ring/gap structure around r s h a dow (Figure[5]), while other models remain largely smooth 
and featureless throughout in both the original and 1/r 2 scaled images. 

Remarks on the implications of the models and observations: 

1. A faint disk in NIR direct imaging observations with a relatively steep surface bright¬ 
ness radial profile (i.e. r~ -2 ) is consistent with small grains being heavily depleted 
and the entire disk is self-shadowed by the inner rim. The surface brightness thresh¬ 
old depends on the stellar luminosity. For a 10 4 K, 2.5 M 0 , and 2 R Q central star, 
the threshold is about l/(r/140 AU) 2 mJy arcsec -2 . Such kind of faint disks may 
have already been found in observations (e.g. HD 139614, Subaru direct imaging, M. 
Fukagawa, private communication, Hashimoto et al. in prep.). 

2. Even with the sharpest possible rims as the parametrized rims in B models, the shadow 
edge fails to manifest itself as prominent rings/arcs in scattered light images. The 
highest ring/gap contrast ratio in full resolution images is only 1.33, and practically no 
ring/gap structure is found in convolved images with an angular resolution of 0.05". 
Note that the convolved images in this study have the ideal chance to review potential 
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disk structures, as (1) they are convolved by a Gaussian PSF, which smears the intrinsic 
full resolution images less than a realistic PSF that has a similar kernel but more 
extended wings; (2) no observational noises and instrument effects are included, which 
are expected to further weaken the detectability of features in images. 


3. Realistic rims in A models can only produce structures even smoother than B models 
in scattered light images. Radial structures with high contrast recently observed in 


NIR direct imaging surveys, such as the arcs in the HD 163296 system (Garufi et al. 


2014), the gap in the HD 169142 system (Momose et al. 2015), the spiral arms in the 


SAO 206462 system ( 

Muto et al. 

2012 

Garufi et al. 

2013 

and the MWC 758 system 

(Grady et al. 2013; 

Benisty et al. 2015) 

, and various rings and arcs in transitional disks 

(e.g. 

Hashimoto et al. 2012; 

Mayama et al. 2012 

) are unlikely to be the shadow edges of 


puffed up rims. It is more likely that they are created by local density or scale height 
structures (Juhasz et al.||2015), such as gaps (e.g. Zhu et ah||2011; Dodson-Robinson & 


Salyk||2011 ; de Juan Qvelar et al.||2013 ; Pinilla et al.||2015 ; Dong et ah|2014[ ) and spiral 


arms (e.g. Dong et al. 2015b) produced by planets in the density distribution^} 

4. Though unable to create arcs/rings and gaps in disks, a shadow cast by a puffed up rim 
may still have an effect on the surface brightness radial distribution of scattered light 
images. With a vertically sharp rim, scattered light images show a 3-stage broken power 
law: 2 steep radial profiles in the inner and other disk connected by a shallow transition 
region around the shadow edge. Such systems may have already been observed: the 
H band scattered light image of TW Hydrae has a shallow r~ 03 transition region 
located at r ~20-40 AU, in between 2 steep r -1 ' 4 and r ~ 2,7 profiles in the inner and 
outer disk (Akiyama et al.|| 2015). The shallow region may actually be a gap, which can 
be either confirmed or ruled out by future observations with better image quality and 
S/N ratio. If the shadow region is not a gap, then comparing with the B models this 
may suggest that the disk is self-shadowed, and the shadow edge is around ~ 30 AU 
(note that TW Hydrae has an r ~ 4 AU inner hole at the center, Calvet et al. 2002 


Hug hes et al. 2007} therefore perhaps the shadow is created by a puffed up cavity wall 
at ~ 4 AU instead of a rim at ~ 0.5 AU). However a self-consistent structure puffed 
by the HSEQ condition is unlikely to account for the observation, as the shadowed- 
to-flared-disk transition is too smooth. The needed puffed up structure has to have 
a somewhat sharp cutoff, perhaps sculpted by other mechanisms in the system. The 
realization of such a structure with a physical origin can be the subject of future studies. 


4 An inclined inner disk can produce azimuthally confined shadows with high contrast on the outer disk, 
as shown in|Marino et al.| (2015). 
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Table 1: Model Properties 


Model 

Parametrized/Natural Rim a 

/grains,small 

-Hrim,opaque ( ) 

^shadow (AU) 

Ring/Gap Contrast 6 

Base 

Natural 

1 

N/A 

N/A 

N/A 

Al 

Natural 

10 -1 

N/A 

N/A 

N/A 

A2 

Natural 

10" 2 

N/A 

N/A 

N/A 

A3 

Natural 

10" 3 

N/A 

N/A 

N/A 

A4 

Natural 

icr 4 

N/A 

N/A 

N/A 

B1 

Parametrized 

1 

10 

5 

N/A 

B2 

Parametrized 

1 

12.5 

13 

1.33/1.03 

B3 

Parametrized 

1 

15 

30 

1.10/1.02 

B4 

Parametrized 

1 

17.5 

52 

N/A 

B5 

Parametrized 

1 

22.5 

147 

N/A 

B6 

Parametrized 

1 

27.5 

N/A 

N/A 


Note. — Properties of the models, (a) Whether the inner rim is a naturally puffed up rim (i.e. the vertical 
structure of the rim is self-consistently determined by the HSEQ condition) or a parametrized artificial opaque 
rim (i.e. the rim structure is fixed and NOT evolved under the HSEQ condition; see Section [ 2 ] for details), 
(b) The mass fraction of the sub-/im-sized small grains in the model ra gra i nS5Sma n/ra gra i nS5 totai (the rest of the 
grains are up to mm-sized big grains), (c) the height of the artificial rim expressed in the angle extended by 
half the rim to the star, (d) The edge of the shadow cast by the rim, defined as the turning point on the 
r = 1 surface (Figure |2|. A distinct shadow edge does not exist in the A models and B6. (e) The contrast 
of the ring/gap feature around r s h a dow in face-on H band PI image, defined as the ratio of the azimuthally 
averaged surface brightness at the peak of the ring to the bottom of the gap inside. Only model B2 and 
B3 has such a feature (i.e. a region with PI^ r >0 ). The 2 numbers are for full resolution raw images (left) 
and convolved images (right, Gaussian PSF FWHM=0.05 // )* ^shadow in B1 is too small to have a meaningful 
ring/gap definition for current NIR imaging facilities. 
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T/(K) T/(K) 



r (AU) log(r/AU) 

Fig. 1.— 2D density (p) and temperature (T) structure as a function of radius r and polar 
angle 9 in spherical coordinate for model Base (left column) and B2 (right column). The 
blue line in each panel marks the disk surface, defined as where the radial optical depth r 
at 0.4 pm at each 6 angle from the star reaches 1 (the stellar flux peaks around 0.4 pm). 
The puffed up rim in model Base and the parametrized artificial rim in model B2 at r su b 
are labeled. The region under the dashed line in model B2 is the umbra, where the entire 
star is not visible to the disk (above this line the polar region of the star may be visible). 
The vertical dotted line in B2 marks the edge of the shadow at r shadow = 13 AU. The r = 1 
surface traces out the boundary between the hot disk atmosphere and the cool interior; the 
umbra region is sufficiently cold so that the disk collapses to the midplane; while the r = 1 
surface and the umbra boundary further separate out a region that receives a small amount 
of direct stellar radiation from the polar region of the star and is moderately puffed up in 
the vertical direction. 
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Fig. 2.— Radial r = 1 surface in spherical coordinates. A shadowed disk has a flat surface 
(constant 0 with increasing r) while a flared disk has a rising surface (increasing 0 with 
increasing r). In the A series the shadow edge is difficult to define, as the disk transits from 
the self-shadowed region to the flared region in a smooth way, due to a lack of sharp rim edge. 
In the B series the transition point is obvious, marked by the dot at [r s h a dow, ^Aim.opaque] on 
each curve, where the surface comes out the shadow cast by the inner rim and starts to rise. 
Model B6 stays in the shadow of its high rim throughout the disk. This surface transition 
point rises and moves outward as the height of the artificial rim increases, and the disk 
surface always converges to model Base beyond r s h a dow, indicating that the flared part of the 
disk in all models are similar. 
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Fig. 3.— Mid-plane temperature profile for the B series. The shadow edge at r shadow is 
marked by the dots on the curves. There are two bumps on the curves. The first bump 
occurs when the disk comes out of the umbra, and the second bump occurs around the 
shadow edge r s h a dow- The transition of T m j, ; i p ] arie around r s h a dow is smooth, without a sudden 
jump. 
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Fig. 4. 


SED of the models at 45° inclination. 
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Fig. 5.— H band polarized intensity images for B2 and B4, showing full resolution 
raw images (left column in each model) and images convolved by a Gaussian PSF with 
FWHM=0.05" (right column in each model) at face-on (top row) and 45° inclination view¬ 
ing angles (bottom row, lower half is the front side). Models are assumed to be at 140 pc 
away. The dashed circle in face-on images marks r s h a dow- The red dot in all panels marks 
the location of the star. The scattered light from the artificial inner rim is removed. The 
central bright region in B4 is masked out for better illustration of the outer disk. The PSF 
size is indicated at the lower right corner in the convolved images. B4 and raw images of B2 
are shown in logarithmic scale, while convolved images in B2 are in linear scale to highlight 
the shadow edge feature. The ring/gap structure at r s h a dow in B2 is marginally traceable, 
while it is not visible in B4. 
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Fig. 6.— Surface brightness as a function of radius for face-on H band PI images for all 
models. Upper panels show the intrinsic surface brightness, and lower panels show the (1/r) 2 
scaled values. Solid curves are convolved images (Gaussian PSF FWHM= 0.05"; targets at 
140 pc); dashed curves are full resolution images; and dotted lines indicate the trend lines. 
v ,shadow is marked for B1 to B5 as a dot. The A models always have a smooth radial profile 
reasonably well characterized by a single power law, while the general surface brightness 
decreases by 1-2 orders of magnitudes from Base to A4.. On the other hand the B series are 
characterized by a 3-stage broken power law, with two steep radial profiles in the inner and 
other disk and a shallow one in between around r s h a dow- Even in the extreme cases with the 
sharpest possible rim edge in the B models, the shadow edges do not manifest themselves as 
prominent rings in images. Merely they create a region in between the inner shadowed disk 
and the outer flared disk with a shallow radial profile close to r°. 












